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Abstract. The outbursts of the dwarf nova WZ Sagittae occur on timescales 
which are much longer than expected from the standard thermal viscous disc 
instability model, which correctly reproduces the outburst behaviour of many 
other dwarf novae. In addition short 'echo' outbursts have sometimes been 
observed after the main outburst. The increased timescales could be explained 
by a magnetic propeller, anchored to the primary star, which clears the inner 
disc. We argue that the propeller, combined with the resonances of the binary 
orbit, may also be the cause of the echo outbursts. We present numerical models 
which reproduce the echo outbursts. 



1. Introduction 

WZ Sagittae is a member of the SU Ursae Majoris class of dwarf novae. These are 
short period systems which exhibit long superoutbursts in addition to ordinary 
outbursts. During these superoutbursts periodic modulations of the optical light 
curve, known as superhumps are observed. These features have a frequency 
which is slightly longer than that of the binary orbit, and which is usually 
attributed to a precessing disc. Superoutbursts in the SU UMa systems typically 
occur with recurrence times of weeks or months. 

Stars from the subclass to which WZ Sge belongs are sometimes known 
as TOADs or Tremendous Outburst Amplitude Dwarf novae. These systems 
have very long recurrence times and correspondingly bright superoutbursts. WZ 
Sge itself showed a 33 year recurrence time, until this pattern was broken by 
the 2001 outburst which occurred after only 23 years. WZ Sge, the brightest 
of all dwarf no v ae, is also notable for its extreme mass ratio. According to 
ISkidmore et"afl (I2002T ) M 2 = 0.045 ± 0.003 and Mi = -74 ± 0.07 M . 
It has a binary period of P Q rb = 81.6 min ( Krzeminskil 1962) a nd a primary 
spin period of P S pin = 27.87 s. In this outburst fe.g iPatterson et al.l (2002) 
superhumps were observed, but pri or to these a m o dulat ion at the frequency 
of the binary orbit was observed bv IPatterson et all (|2002T ) as discussed by e.g. 



iKunze fc Sneithl (|2004ft . In addition, the main outburst was followed by twelve 



echo outbursts. Echoes have also been ob served in the similar system EG Cancri 
( Osaki. Mever Mever-Hofmeisterll200ll h 
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Numerical mode l s inco rporating the thermal viscous disc instability model 
(DIM) due to iHoshil ( 19791 ) . have reproduced the outburst behaviour of many 
dwarf noave. These models have not however, without using extremely low 
viscos ities, been ab le to explain the long recurrence time of WZ Sge type sys- 
tems. lO saki (1989) extended the DIM to binaries with small mass ratios when 
IWhitehurstl (jl988T ) showed that discs in such systems could become tidally unsta- 
ble. Osaki proposed that, during the normal outburst cycle, material is allowed 
to build up in the outer disc. When the disc becomes sufficiently large that 
it may access the 3:1 resonance of the binary orbit, it will become eccentric, 
enhancing dissipation a nd the loss of angular momentu m, leading to a superout- 
burst. Simulations (e.g. iTruss. Murray k, Wvnnll200lh supported this theory. 

The extended recurrence times of the WZ Sge type systems may be ex- 
plained by the use of a magnetic propeller to evacu ate the inner regions of the 
disc. T his mechanism, similar to that proposed bvfMatthe ws. Speith &; Wvnnl 
2004a|) for the FU Orionis stars, is detailed in the forthcoming paper bv lWvnn et aD 
The magnetic propeller is a mechanism by which the inner regions of 



an accretion disc can be evacuated. The co-rotation radius R co is the radius at 
which a Keplerian orbit has the same angular frequency as that of the star's spin. 
If the primary magnetic field rotates as a solid body then this radius can be said 
to divide the magnetic accretion and magnetic propeller regimes. Particles at 
radii smaller than the co-rotation radius lose angular momentum as a result of 
their interactions with the stellar magnetic field, and fall rapidly onto the star. 
This magnetically enhanced accretion occurs because the particles are moving 
more rapidly than the magnetic field. Elements of plasma which exist outside co- 
rotation are, conversely, accelerated by the magnetic field, and gain angular mo- 
mentum at the star's expense. This is the magnetic prop eller regime. An exam- 
ple o f a strong magnetic propeller is that of AE Aquarii jWvnn, King Sz Hornd 
in which accretion is completely prevented by the magnetic field. What is 
proposed here is a far weaker propeller, driven by a rapidly spinning, moderately 
magnetic white dwarf. Simulations us in g weaker fields have been shown to warp 
accretion discs (|Q' Sullivan et alJl2004T l. lLasota. Kuulkers fc Charlesl (|l999h also 

used a magnetically truncated disc in their model of WZ Sge. 

The cause of the echo outbursts is as yet uncertain. lOsaki. Mever k, Mever-Hofmeisten 

(2001) suggested that the mass which accumulates near the 3:1 resonance in the 
disc is released gradually after the main outburst, leading to echo outbursts. In 
this paper we model a full outburst of WZ Sge using smoothed particle hydro- 
dynamics (SPH). The rapidly spinning white dwarf creates a magnetic propeller 
which, in addition to evacuating the inner disc, aids the creation of a reservoir 
of material in the outer disc. Echo outbursts arise naturally from this configu- 
ration. 



2. Numerical Method 

The simulations in this paper were performed using a fully three dimensional 
SPH code with viscosity switching and a magnetic acce l eratio n as de scribed be- 
low. T he specific code is fully described bv lTruss et al.l (|2000h whilst iMona ghan 
(1992) provides a general review of SPH. SPH is a Lagrangian technique for the 
modelling of fluid dynamics. It works by discretising the fluid into a series of 
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particles, normally with equal mass, which represent fluid elements. These par- 
ticles carry values for fluid properties such as density. The properties of the fluid 
may then be interpolated for any point in space by taking a weighted average of 
the properties of nearby particles. This is done using a smoothing kernel, which 
usually takes the form of a cubic spline function. 

The viscosity in the disc is an artificial viscosity which is introduced into 
the velocity and energy equations of the SPH formulation. The extra term is 
composed of shear and bulk viscosity elements. However ^provi ding that the local 
velocity divergence is negligible, it c an be shown jMurraylll996t) that the viscosity 
acts in a manner very similar to the lShakura &: Sunvaevl ( 1973) a shear viscosity. 



The SPH viscosity parameter £, which is proportional to a, is switched between 
two values, mimicking the behaviour of the DIM. Switching occurs at two trigger 
surface densities. Any annulus of the disc which exceeds Sj^gk will begin to have 
its £ gradually altered to the high state value. Conversely, an annulus which falls 
below Sgw wm have its C value slowly reduced to the lower, quiescent value. 
This is discussed bv iTruss et al.l ( 2000fi . The p ositions of the triggers, which 



are ro ughly proportional to radius, are dervied bv lCannizzo. Shatter fc Wheeler! 
1 19881) . 



Dendv (|l990l ) expresses the magnetic tension acceleration as 

(1) 



1 ( BzB^ 
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where p is density and r c represents the radius of curvature of the magnetic 
field lines. The azimuthal and vertical components of the magnetic field are 
represented by Bj, and B z respectively. We will make the approximation that the 
radius of cu rvature of the field lines is of t he order of the disc scale-height thus, 
r c ~ H (e.g. iPearson. Wvnn Kindll997l ). The ratio of vertical and azimuthal 
field strengths is related to the shear between the disc and the magnetic field. 
If it is assumed that the fie ld rotates rigidly with the star then this ratio can be 
expressed in the form (e.g. lLivio fc Pringlelll9 92) 

B <t> ^ (f*k ~ fk) , 2 s 
B z S7 k 

where £1* indicates the angular frequency of the accreting star while f2 k repre- 
sents the Keplerian angular frequency at a given angle. Equations (^Q) and © 
are combined. Additionally it is assumed that the field is dipolar so that the 
magnetic field can be related to the magnetic moment and distance from the 
star by the expression B z ~ \fj,\ R~ 3 where i?* represents the radius of the star. 
It is also assumed that surface density can be approximated by E ~ pH, where 
H is the scale height of the accretion disc. This allows us to make the following 
approximation 

Bl (n k -n*)_ ^i^K-p kmR -6 (^- v *) (3) 



6 AnpH n k 4vrS u k u k 

In this foundation the acceleration falls to zero when i> k = v+ at co-rotation as 
expected. It may be useful to compare the value constant k m with the (3 used 
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Figure 1. Upper frame: Light curve of WZ Sge in 2001, showing a rapid 
rise and subsequent decay punctuated by an episode of 12 echo outbursts. A 
freehand cu rve has been added to improve visibility. This article originally 
appeared in lPatterson et aP l)2002|) . Copyright 2002, Astronomical Society of 
the Pacific; reproduced with permission of the Editors. Lower frame: Results 
of simulations showing the occurrence of echo outbursts. 



in iMatthews et aT | (|2004bT) . It can be shown that 

where G is the universal constant of gravitation and M* is the mass of the star. 



3. Results 



Simulations were performed using the code described above. The results pre- 
sented used a very high viscosity, such that the Shakura-Sunyaev alpha parame- 
ters in the hot and cold states are 0.6 and 0.06 respectively. The triggers are also 
placed rather close together and these initial calculations were performed using 
a low resolution of 30,000 particles, since they were part of a large parameter 
study. The very high viscosity has the effect of speeding up the viscous processes 
and hence reducing processing time. The viscosity switching triggers were acti- 
vated only after the disc had been filled. A weak magnetic field of B ~ 1 x 10 5 G 
was used. As shown in Figure 1, the main outburst in our simulation is followed 
by a series of echo outbursts which are of a similar form to those in the obser- 
vation which is also plotted. The timescale for the echo outbursts is an order of 
magnitude shorter than that observed. This is as expected since physical disc 
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Figure 2. Upper frames: Snapshots of the accretion disc from our simula- 
tions showing log surface density and the viscous parameter £ — 10a respec- 
tively. The binary separation is represented by a in both plots. Lower frame: 
A plot of azimuthally averaged surface density as a function of radius also 
showing the viscosity switching triggers. 



alphas are thought to be of order 0.1 and 0.01. The echoes do not at present 
switch off, but continue up to the end of the run. 

Figure 2 shows the state of the disc during the echo outburst phase. The 
upper panels show surface density and the viscous alpha in the disc. The two 
circles represent the 2:1 and 3:1 orbital resonances. The lower panel is a plot 
of the azimuthally averaged surface density as a function of radius. The central 
hole in the accretion disc, due to the magnetic propeller, is visible in all these 
plots. This latter plot also illustrates the positions of the two viscosity triggers 
in the simulation. Together the plots show how the echo outbursts occur in the 
simulations. After a main outburst, some material remains in the outer disc. 
Spiral waves are excited, as a result of the binary potential. The associated 
density enhancements, when augmented by continuing mass transfer, aid the 
transport of matter towards the centre where, as accretion is inhibited by the 
magnetic propeller, an echo outburst occurs. The echo fades when the propeller 
is overcome by the viscous forces and accretion onto the star resumes. Echoes 
do not appear if the simulation is performed without the magnetic field. 
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4. Discussion 

In this paper echo outbursts are reproduced in an SPH model of WZ Sge. These 
echoes are similar to those observed in the light curve of the system However, 
it is not yet clear whether this mechanism is likely to be the one which occurs 
in reality. In order to establish this there are a number of refinements which 
will be made to the model. Firstly the viscous alphas used in the model will be 
brought closer to their accepted values. It is not desirable to rely on time scaling 
arguments in this system since neither magnetic nor resonant effects scale with 
alpha. The ~ 30, 000 particles used are not sufficient to fully resolve such a disc 
in three dimensions, so higher resolution calculations will be made. In addition 
a mechanism for the termination of the echoes will be sought. 
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